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A-layer transition metal dichalcogenides provide a unique platform to investigate the evolution of 
the physical properties between the bulk (three dimensional) and monolayer (quasi two-dimensional) 
limits. Here, using high-resolution micro-Raman spectroscopy, we report a unified experimental 
description of the F-point optical phonons in A-layer 2H-molybdenum ditelluride (MoTe 2 ). We 
observe a series of A-dependent low-frequency interlayer shear and breathing modes (below 40 cm“^, 
denoted LSM and LBM) and well-defined Davydov splittings of the mid-frequency modes (in the 
range 100 — 200 cm“^, denoted iX and oX), which solely involve displacements of the chalcogen 
atoms. In contrast, the high-frequency modes (in the range 200 — 300 cm“^, denoted iMX and 
oMX), arising from displacements of both the metal and chalcogen atoms, exhibit considerably 
reduced splittings. The manifold of phonon modes associated with the in-plane and out-of-plane 
displacements are quantitatively described by a force constant model, including interactions up to 
the second nearest neighbor and surface effects as fitting parameters. The splittings for the iX and 
oX modes observed in A-layer crystals are directly correlated to the corresponding bulk Davydov 
splittings between the E^njEig and BiujAig modes, respectively, and provide a measurement of 
the frequencies of the bulk silent E 2 u and optical phonon modes. Our analysis could readily be 
generalized to other layered crystals. 

Keywords: Two-dimensional materials, layered crystals, transition metal dichalcogenides, 
MoTe 2 , Raman spectroscopy, interlayer breathing and shear modes, force constants, Davydov split¬ 
ting, surface effects. 


Introduction In the wake of graphene, a vast fam¬ 
ily of layered materials is attracting tremendous atten¬ 
tion [T]. Now available in the form of A-layer crystals, 
the latter exhibit peculiar physical properties that com¬ 
plement the assets of graphene and offer exciting perspec¬ 
tives to design van der Waals heterostructures [T]. Semi¬ 
conducting transition metal dichalcogenides (MX 2 , with 
M = Mo, W and X = S, Se, Te) are among the most ac¬ 
tively investigated layered crystals [5]. Indeed, although 
bulk MX 2 exhibit indirect bandgaps, monolayer MX 2 are 
direct bandgap semiconductors [a n with remarkable 
spin, valley [5] and optoelectronic properties [6]- More 
Generally, A-layer MX 2 crystals provide an ideal plat¬ 
form to uncover the impact of symmetry breaking and 
interlayer interactions on the electronic, optical and vi¬ 
brational properties, from the bulk (three-dimensional) 
to the monolayer (quasi two-dimensional) limit. 

In particular, in A-layer MX 2 , interlayer interac¬ 
tions result in a splitting of all the monolayer phonon 
modes [THIS] (see Table [B. The latter effect is known 
as the Davydov splitting |I7| and is closely related to 
the force constants that govern the vibrational properties 
of MX 2 |9]. The Davydov splitting has been previously 
studied in polyaromatic molecules [15], thin films m, 
and bulk layered crystals, including MX 2 [THSIIISIIII]- 
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Recently, low-frequency Raman spectroscopy has been 
employed in A-layer M 0 S 2 [HI and WSe 2 [M] 

to uncover the fan diagrams associated with the frequen¬ 
cies of the interlayer shear (LSM) and breathing (LBM) 
modes, which arise from the Davydov splitting of the zero 
frequency acoustic modes (see Fig. [^d)). Related split¬ 
tings of the higher frequency modes involving intralayer 
atomic displacements have been reported in A < 5-layer 
MoSe 2 [HI HT] and WS 2 [Hj- Thus far, such splittings 
have solely been reported for the out-of-plane Raman- 
active phonon with Aig symmetry in the bulk, whereas 
other high-frequency modes, such as the in-plane phonon 
with E 2 g symmetry in the bulk, exhibit anomalous A- 
dependent frequency shifts [8HTU1 [HI IHl [HH?T] but no 
splitting. However, a unified description of the Davydov 
splitting in A-layer MX 2 is still lacking. 

In this letter, we quantitatively investigate the Davy¬ 
dov splitting of all the F-point optical phonon modes, 
over the range 4 — 300 cm“^ in A-layer 2H molybdenum 
ditelluride (MoTe 2 ) crystals. As in other MX 2 , A-layer 
MoTe 2 displays a set of [A/2J (herein [J ([]) denote the 
floor (ceil) functions that map A/2 to the largest previ¬ 
ous (smallest following) integer) well-separated LSM and 
LBM below 40 cm“^. In the range 100—200 cm“^, the in¬ 
plane and out-of-plane modes involving solely displace¬ 
ments of the chalcogen (Te) atoms (herein denoted iX and 
oX) exhibit pronounced splittings into L-^/2j and \N/2] 
Raman-active modes, respectively. The Davydov split¬ 
ting between the bulk in-plane Eig Raman-active and 
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FIG. 1. (a) Side and top view of the crystal structure of 2iF-MoTe2. The red rhombus represents the unit cell, (b) Optical 
image of a A^-layer MoTe 2 crystal deposited onto a Si/Si02 substrate, (c) Raman spectra of monolayer, bilayer and bulk 
MoTez in the parallel (XX, thick colored solid lines) and perpendicular (XY, thin grey lines) polarization configuration. The 
spectra are vertically offset for clarity and the asterisk highlights residual contributions from the exciting laser beam, (c) 
Atomic displacements and irreducible representations associated with the F point phonon modes in monolayer, bilayer and 
bulk MoTez. The Raman (R) and/or infrared (IR) activity are indicated, and stars denote silent modes. The zero frequency 
acoustic (A) modes (LA, TA, ZA) and their irreducible representations are also shown for clarity. 


E 2 u silent modes (out-of-plane Aig Raman-active and 
Biu silent modes) is 2.7 cm“^ (4.7 cm“^). In the range 
200 — 300 cm“^, the modes involve motion of both the 
metal (Mo) and chalcogen (Te) atoms. The in-plane iMX 
mode exhibits no observable splitting. The out-of-plane 
oMX mode exhibits a small splitting of approximately 
1 cm“^, which is actually not a Davydov splitting but a 
splitting due to surface effects. 

All the observed trends are fitted using a finite linear 
chain model that includes inter- and intralayer force con¬ 
stants up to the second nearest neighbor, as well as sur¬ 
face effects at the edges of the chain . This model nat¬ 
urally explains why the interlayer force constants, which 
directly determine the manifold of rigid layer modes, also 
give rise to sizable Davydov splittings for the iX and oX 
modes but to reduced splittings for the iMX and oMX 
modes. Conversely, surface effects are responsible for the 
apparent downshift of the iMX mode as N augments and 
for the slight splitting of the oMX mode. 


In its stable form, bulk MoTez displays a trigonal pris¬ 
matic {2H polytype) structure with AbA/BaB stacking 
(see Fig. and shares many properties with the widely 
studied M 0 S 2 , MoSez, WS 2 and WSez crystals [H IT^ . 
MoTez has recently attracted particular interest due to 
its lower transport bandgap [33U36] and near-infrared 
emission [37113H]- Indeed, at room temperature, mono- 
layer MoTez exhibits a direct optical bandgap at 1.1 
eV, whereas the bulk material has an indirect optical 
bandgap slightly below 1.0 eV [2l[37l[38]. Bulk MoTez 
exhibits 3 acoustic and 15 optical phonon modes at F as 
is the case for MoSz and MoSez [Ml [IMS]. As illus¬ 
trated in Fig[^ the manifold of optical phonons breaks 
down into 5 (doubly degenerate) in-plane modes and 
5 out-of-plane modes. According to group theory [53] . 
the integrated intensity of the in-plane Raman-active 
modes is expected to be independent on the angle 9 be¬ 
tween the polarization of the incoming and scattered pho¬ 
tons, whereas the integrated intensity of the out-of-plane 
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TABLE I. Irreducible representations of the optical phonon modes at T for single-, bi-, N- layer and bulk MoTe 2 . Bold characters 
denote Raman-active modes in a backscattering geometry. Note that modes with Eig and E" symmetry are Raman-active but 
not observable in a backscattering geometry m and that modes with E' are Raman- and infrared-active. Stars (*) denote 
silent modes. All the other modes are infrared active. 


modes should exhibit a cos^ 9 dependence. 

Results Figure [2 shows the micro-Raman spectra of 
monolayer, bilayer and bulk MoTe 2 recorded at Ei, = 
2.33 eV for 0 = 0 (parallel configuration, XX) and 
6 = TTj2 (cross-polarized configuration, XY). We first 
address the low-frequency range below 40 cm“^. In bulk 
MoTe 2 , we observe only a single feature at 26 cm“^ that 
shows similar intensities in the XX and XY configura¬ 
tions. The latter is assigned to the in-plane interlayer 
shear mode (LSM) with E 2 g symmetry [5D]. For N = 2, 
we observe a prominent feature at 28 cm“^ that shows 
strong extinction in the XY polarization and a fainter 
feature at 18 cm“^ whose intensity is similar in the XX 
and XY configurations. The former is thus assigned to 
the out-of-plane layer breathing mode (LBM) with Aig 
symmetry, whereas the latter is assigned to the LSM with 
Eg symmetry. The LBM has B 2 g symmetry in the bulk 
and is silent. As expected, we do not observe any inter¬ 
layer mode for Af = 1. 

In the mid- (100 — 200 cm“^) and high-frequency 
(200 — 300 cm“^) ranges, the Raman spectra of bilayer 
MoTe 2 displays four one-phonon features, which have 
previously been identified as originating from the fol¬ 
lowing intralayer displacements: (i) the in-plane, out- 
of-phase vibration of the Te planes, with Eig symmetry 
in the bulk (iX mode, near 120 cm“^), (ii) the out-of- 
plane, out-of-phase vibration of the Te planes, with Aig 
symmetry in the bulk (oX mode, near 170 cm“^), (iii) 
the in-plane vibration of the Mo and Te planes against 
each other, with E 2 g symmetry in the bulk (iMX mode, 
near 230 cm“^), and (iv) the out-of-plane vibration of the 
Mo and Te planes against each other, with B 2 g symme¬ 


try in the bulk, (oMX mode, near 290 cm“^) [13 UHl ES] ■ 
The bulk iX and oMX modes are predicted to be Raman 
inactive in a backscattering geometry and silent, respec¬ 
tively. However, both modes appear as faint features in 
thick MoTe 2 flakes (considered as a bulk reference). This 
surprising observation, also reported recently on other 
MX 2 might be a consequence of the finite penetration 
depth of our laser due to the strong optical absorption of 
MoTe 2 0 E7] or may arise from a breakdown of the Ra¬ 
man selection rules due to resonance effects [13 IHl |40] . 
As predicted by group theory [T^HTl] . the iX and oMX 
modes are not observed in monolayer MoTe 2 (see Figs. 
andj^a)). We verified that the oX and oMX features 
nearly vanish in the XY configuration, whereas the in¬ 
tegrated intensities of the iX and iMX features do not 
change (see Fig. [^. 

Figure shows the evolution of the low-frequency 
modes (LSM and LBM) from N = 2 to N = 12- 
layer MoTe 2 . As previously reported in multilayer 
graphene gll HU and recently in other MX 2 , [M m- 
EHIIZI, a set of Y-dependent low-frequency Stokes and 
anti-Stokes peaks appears for Y > 2. The number of 
detected peaks increases with Y and the peaks can be 
separated into branches that seemingly stiffen or soften 
with increasing Y (see the dashed lines in Figj^a)). In¬ 
terestingly, compared to a reference recorded in the XX 
configuration, the integrated intensity of peaks belonging 
to a branch that softens with increasing Y drops by more 
than one order of magnitude in the XY configuration, 
whereas the integrated intensity of the peaks that belong 
to a branch that stiffens with increasing Y is marginally 
affected. Therefore, the branches that soften (stiffen) 
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FIG. 2. (a) Polarization-resolved low-frequency Raman spectra of = 1 to = 12 layer MoTe 2 and of bulk 2Jf-MoTe2, 
in the parallel (XX) and transverse (XY) configuration recorded at Sl = 2.33 eV. The interlayer breathing modes (LBM) 
largely dominate the spectra in the XX configuration and their intensity is reduced by more than one order of magnitude in 
the XY configuration. The shear modes (LSM) are not sensitive to the polarization configuration. The dashed lines follow 
the frequencies of each LSM and LBM (see also Fig. [^. (c-d) Polarization-resolved low-frequency micro-Raman spectra of 
(b) N = 9 and (c) Y = 10 layer MoTe 2 in the parallel (XX) and transverse (XY) configuration. The four (five) expected 
Raman-active LBM (denoted Ba to Bd (Ba to Be) for Y = 9 (Y = 10) and two LSM (denoted Sa and Sb) are observed. The 
asterisks highlight residual contributions from the exciting laser beam. 


with increasing N are assigned to the LBM (LSM). We 
are able to resolve two branches of LSM and, remark¬ 
ably, five branches of LBM, i.e., the complete manifold 
of Raman-active LBM up to Y = 11 (see Table [l] and 
Fig. [5]). 

The Raman spectra of the mid- (iX, oX) and high- 
frequency (iMX and oMX) modes in Y-layer MoTe 2 
are shown in Fig. For Y > 3-layer MoTe 2 , we ob¬ 
serve a prominent splitting of the oX-mode feature (see 
Fig.^b)), whereas in the bulk limit, one recovers a sin¬ 
gle symmetric feature (assigned to the Raman-active Aig 
mode). Interestingly, as shown in Fig. [^a), the iX-mode 
feature also splits, but only for Y > 4. We can resolve up 
to three subfeatures for Y = 6 and 7, but the Raman sig¬ 
nal in Y > 8-layer MoTe 2 becomes too small to perform 
a quantitative analysis. At Yl = 1.96 eV, the oMX-mode 
feature also exhibits a modest splitting, on the order of 
1 cm“^, for Y > 4. Two subpeaks can be distinguished. 
However, the evolution of their frequencies does not fol¬ 
low a specific trend as a function of Y (see Fig. id) and 
Fig. I^f)). In contrast, the iMX-mode feature exhibits a 
faint shoulder on its high-energy side (see Fig. ic)), but 


no appreciable splitting can be resolved. However, the 
iMX feature downshifts as Y increases, consistently with 
previous reports on Y-layer MX 2 [HHlQl CSl ESI 1^ Ei- 

1311 ISZl- 

The splitting of a bulk phonon mode in a Y-layer 2H- 
MX 2 can be understood based on a group theory anal¬ 
ysis [SHISl HE]- Bulk MoTe 2 belongs to the Deh non- 
symmorphic space group and crystals with odd or even 

Y belong to the D^h and D^d symmorphic space groups, 
respectively. In bulk 2 Y-MX 2 , the (silent) LBM and (Ra¬ 
man active) LSM correspond to an out-of phase displace¬ 
ment of adjacent layer pairs, whereas the iX, iMX, oX, 
and oMX atomic displacements give rise to two modes, in 
which the relative motion of equivalent atoms belonging 
to two adjacent layers is either in-phase or out-of-phase 
(see Fig. Bd))- In a Y-layer system, interlayer inter¬ 
actions and symmetry breaking result in a manifold of 

Y - 1 LSM and LBM, and Y iX, oX, iMX and oMX 
modes. Table |T] summarizes the irreducible representa¬ 
tions and the activity of these optical phonon modes for 
mono-, bi-, Y-layer, and bulk MoTe 2 . 

In agreement with group theory predictions, our Ra- 
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FIG. 3. Normalized micro-Raman spectra of the (a) iX, (b) oX, (c) iMX, and (d) oMX mode-features in N-layer MoTe 2 . The 
laser phonon energy i?L and the elementary intralayer displacements are indicated. The measured Raman features (symbols) 
are fit to Voigt profiles (solid lines). For the modes that show a Davydov splitting, each subpeak is represented with a colored 
dashed line. A featureless background (grey dashed line) has been considered when necessary. 


man measurements reveal that (i) secondary LSM and 
LBM features appear from N > A (ii) the oX feature 
splits for > 3 and we observe \N/2\ Raman active 
subfeatures in iV-layer MoTe 2 , (iii) the splitting of the 
iX feature into [1V/2J Raman active subfeatures is ob¬ 
served from > 4, (iv) a splitting of the oMX feature 
appears for Ai > 4. The frequencies of all the measured 
Raman features, extracted from Voigt fits (see Figs. 
are reported in Fig. 


Force constant model We now introduce a semiem- 
pirical model to rationalize the observed splittings. Al¬ 
layer MoTe 2 is modeled as a finite linear chain composed 
of N Mo atoms and 2N Te atoms. Within one MoTe 2 
unit nearest neighbor Mo and Te atoms and the pair 
of second nearest neighbor Te atoms are connected by 
springs with force constants per unit area a and S respec¬ 
tively [S] (see Fig. I^a)). Interlayer interactions are then 
described by two force constants per unit area, /3 and 7 , 
between nearest neighbor Te atoms belonging to adja¬ 
cent layers and between second nearest neighbor Mo and 
Te atoms, respectively. In addition, finite size effects are 
known to lead to a slight reduction of the metal-chalcogen 


bond length on the outer layers (“surface effects”) [3] . As 
a result, to improve our model, effective values Oe > a 
and Se > S are phenomenologically considered for the first 
and A^**' layer. We note that our choice of using the same 
value of Qfe at both ends of the chain is consistent with the 
fact that no significant substrate-induced frequency shifts 
have been observed on the Raman response of MoTe 2 and 
other MX 2 [21 dni HSl Uni [29] . The normal modes = 
(“1,1- “ 2 , 1 - “3,1- ■■■’ “ 2 ,]V- “ 3 ,iv), with 

the normal displacement of the z*** atom (z = 1, 3 for Te 
and z = 2 for Mo) in the j**' MoTe 2 layer (j £ |1,A^]) 
associated with the normal mode k £ |1,3A^] (see Fig. 
I^a)), and the eigenfrequencies of a A^-layer system 
can be obtained from Newton’s equations of motion, 
which lead to the secular equation 

u}lU^=VU'^, ( 1 ) 

involving the 3N x 3N dynamical matrix T) (see Support¬ 
ing Information). 

For in-plane (out-of-plane) displacements, Eq. 0 pre¬ 
dicts three manifolds of N normal modes that correspond 
to (i) the low frequency LSM (LBM) (including the zero 
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FIG. 4. (a) Schematic of the finite linear chain model discussed in the text. Calculated normal displacements associated with 
all the in-plane (b) and out-of-plane (c) optical phonon modes in = 4 and = 5 layer MoTe 2 . The size of the arrows is 
proportional to the amplitude Uij of the normal displacement obtained from the solution of Eq. 0. At a given N, the mode 
frequencies increase from left to right. The irreducible representation of each normal mode is indicated (see also Supporting 
Information). 


frequency acoustic mode) (ii) the mid-frequency iX (oX) 
modes and (iii) the high-frequency iMX (oMX) modes. 
Figure 0b) represents the calculated iX and oX normal 
modes ^sing the parameters in Table 0 and their irre¬ 
ducible representations for A = 4 and N = 5. For the 
low- (LSM, LBM) and mid-frequency (iX, oX) modes, 
the eigenfrequencies increase (decrease) as the layers ex¬ 
hibit more out-of phase (in-phase) relative motion, up to 
the limit of the highest- (lowest-) frequency mode, which 
corresponds to an out-of phase (in-phase) oscillation for 
all layers. This trend is as expected from classical theo¬ 
ries of coupled oscillators. In particular, we can readily 
conclude that the dominant LBM feature corresponds, 
for even N to the out-of-phase oscillation of two blocks 
composed of N/2 layers that vibrate in-phase; for odd N 
to the out-of-phase oscillations of two blocks composed of 
(A — l)/2 layers that vibrate in-phase, whereas the cen¬ 


tral layer stays at rest. In contrast, the dominant LSM 
feature corresponds to an out-of-phase displacement of 
adjacent layers as in bulk crystals nasa (see also Sup¬ 
porting Information). 

As shown in Fig. 0 the fan diagrams associated with 
the in-plane LSM, iX- and iMX-mode frequencies on the 
one hand, and out-of-plane LBM, oX- and oMX-mode 
frequencies on the other hand are very well described by 
the force-constant model. The force constants (per unit 
area) used as fitting parameters and the corresponding 
bulk frequencies (see Supporting Information for their 
analytical expressions) are reported in Table [n| We find 
that all the force constants correspond to a restoring 
force, except for the in-plane second nearest neighbor 
force constant between Te pairs. The values obtained 
here are qualitatively similar to the values that Luo et 
al. reported for M 0 S 2 by fitting the bulk frequencies 
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(obtained from calculations based on the local density 
approximation (LDA)) to a force constant model [9]. 



in-plane 

out-of-plane 

q( 10^® N/m®) 

105 

159 

ae(10^® N/m®) 

107 

163 

d(10i® N/m®) 

2.28 

5.61 

7 ( 10 ^® N/m®) 

0.585 

1.11 

5(10^® N/m®) 

-4.53 

19.8 

(5e(10i® N/m®) 

-4.22 

20.4 


26.8 (A 29 ) 

40.3 (B 29 ) 

^mid(cm"^) 

117.2 {E2 u) 

168.9 (Ri„) 

^iid(cm"^) 

119.9 {Elg) 

173.6 (A 19 ) 


TABLE II. Force constants per unit area and corresponding 
bulk frequencies of the low-frequency (LSM, LBM) and mid¬ 
frequency (iX and oX) modes extracted from the fit of our 
experimental data to the finite linear chain model (see Eq. Q 
and dashed lines in Fig. [^. The irreducible representations 
of the bulk phonon modes are indicated. 

Discussion First, we note that the low-frequency 
branches of LSM and LBM can also be separately mod¬ 
eled using a linear chain of Al-oscillators. Indeed, in first 
approximation, van der Waals interactions between ad¬ 
jacent layers are sufficient to accurately describe the se¬ 
ries of interlayer modes without further consideration of 
the in-plane crystal structure and intralayer force con¬ 
stants [m El EH EU [43]. The LSM and LBM frequen¬ 
cies are then very well approximated by Wiow.fe (N) = 

— cos , with k e |2,AI] (k = 1 cor¬ 
responds to the acoustic mode at = 0). The 

observed Raman-active modes correspond to branches 
with {i) k = N, N — 2 for the LSM Sa and Sb and 
(ii) k = 2, 4, 6, 8, 10 for the LBM Ba to Be (see 
Fig. i (b),(c)). Using the complete model, the bulk fre¬ 
quencies with ^ the MoTe 2 mass per 

unit area, then allow to determine effective interlayer 
force constants per unit area of /3- 1-27 = 3.5 x 10^® N/m^ 
and j5 + 2'y = 7.8 x 10^® N/m^ for the LSM and LBM, re¬ 
spectively (see Table 0. These values are close to those 
derived from the LBM and LSM in M 0 S 2 [H El E5] 
and WSe 2 m- We note that surface effects only affect 
intralayer force constants and have therefore a negligible 
influence on the rigid layer modes. 

Second, a force constant model restricted to the first 
nearest neighbor interactions {i.e., j = S = 0) suffices to 


fit the frequencies of the low- and mid-frequency modes, 
but would then fail to predict the frequencies of the oMX 
and iMX modes. Indeed, Table [IT] reveals that the second 
nearest neighbor interlayer ( 7 ) and intralayer (J) force 
constants are of the same order of magnitude and larger 
than the nearest neighbor interlayer force constant /3, re¬ 
spectively. Using the complete model, the bulk Davydov 
splittings of the iX and oX modes can conveniently be 

expressed as w+^ ^F Mx 

the chalcogen (Te) mass per unit area, leading to the val¬ 
ues of 2.7 cm“^ and 4.7 cm“^ for the iX and oX modes, 
respectively (see Supporting Information). 

According to our calculations and to group theory m 
[in EH] (see Fig. |1|4| and Table |T|, the mid-frequency 
Raman-active modes that are observable in our backscat- 
tering experiments correspond (i) to the second-lowest, 
fourth-lowest,... frequency mode for the iX phonons and 
(ii) to the highest, third-highest... frequency mode for 
the oX phonons. These distinct symmetry properties re¬ 
sult in a set of softening and stiffening branches in the 
experimentally measured fan diagrams in Fig. [^c) and 
Sd), respectively. A remarkable validation of this sym¬ 
metry analysis is that the highest frequency iX mode that 
can be observed has the highest frequency {Eg symme¬ 
try) for even N and the second highest frequency {E' 
symmetry) for odd N. As a result, the iX-mode fre¬ 
quency for AI = 3 is lower than for N = 2, and the two 
observed iX-mode subfeatures for Al = 4 are slightly up- 
shifted relative to their counterparts recorded for N = 5 
(see Figs. |^and[^. We also note that our model pre¬ 
dicts a somewhat lower frequency for the oX (A() mode 
in monolayer MoTe 2 (see Fig. [^e)). This discrepancy is 
due to a stronger bond length contraction in the limit 
of a monolayer [ 3 , as compared to the bilayer or to the 
outer layers in Al > 2 -layer samples. 

Third, we focus on the high-frequency iMX and oMX 
modes. Our model predicts very small anomalous Davy¬ 
dov splittings (below 0.2 cm“^) for the iMX and oMX 
modes (consistent with previous studies of the iMX mode 
in bulk MX 2 (HHHEnil^) and a critical influence of sur¬ 
face effects. Indeed we find that among the N iMX and 
oMX normal displacements, two modes predominantly 
involve atomic displacement of the first and/or layer 
and can be regarded as surface modes. Such surface 
modes are not predicted for the iX and oX displacements 
(see Fig. and the Supporting Information). As a re¬ 
sult, the surface terms Oe and Se, will only bring a small 
correction to the fan diagrams associated with the mid¬ 
frequency modes, but will shape the manifold of high- 
frequency phonon modes. For N > 3, our model predicts 
a fan diagram for N — 2 quasi-degenerate iMX and oMX 
modes, where atomic displacement chiefly occurs in the 
inner layers and two surface modes (see Fig. and 
Fig. c),(f). The strength of the surface effects will set 
the frequencies of the surface modes relative to that of 
the other modes. 

A tentative explanation for the pronounced downshift 
















N 


N 




N 


N 


FIG. 5. Fre quencies of the (a) LSM (b) LBM (c) iX, (d) oX, (e) iMX, (f) oMX modes extracted from fits of the spectra displayed 
in Figs 2|3 (black circles) as a function of the number of layers N. The red open dries in (a)-(f) are the frequencies of the 
Raman-active modes calculated by solving Eq.Q with the fitting parameters in Table [III The grey-filled triangles in (a), (c), (f) 
(resp. the open squares in (a)-(f)) are the frequencies of the E" modes that are not Raman active in a backscattering geometry 
(resp. of the infrared-active modes) also predicted by Eq. 0 using the same fitting parameters. The solid and dashed lines 
connect the calculated frequencies and are guides to the eye. The upper and lower horizontal dotted lines in (a)-(b), (c)-(d), 
and (e)-(f) correspond to the bulk frequencies and respectively (see also Table fidl. The green open circles 

(squares) in (e) and (f) correspond to the Raman-active (infrared-active) surface modes (see also 0b)). 


observed for the iMX mode is then that the surface mode 
dominates for small N and that the modes arising from 
the inner layers provide most of the Raman intensity as 
N increases. This scenario is consistent with calculations 
based on a bond polarizability model [H] , which was also 
recently used for layered BN and BN nanotubes [IS] and 
which is explained in the Supporting Information. The 
model predicts that the surface mode will dominate in 


the iMX peak for iV < 5. Let us also note that the iMX- 
mode feature in TV-layer MoTe 2 is appreciably broader 
(with a full width at half maximum (FWHM) decreas¬ 
ing from 3.5 cm“^ for TV = 1 down to 2.6 cm“^ in the 
bulk) than the iX, oX and oMX peaks, whose FWHM 
are approximately 1 cm“^. This broadening prevents us 
from unveiling any splitting of the iMX-mode feature in 
iV-layer MoTe 2 . This splitting, according to our force- 
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constant fit is very small anyway (see Fig.[^(e)) and also 
our ab-initio calculations for phonons of bulk MoTe 2 yield 
a minimal anomalous Davydov splitting of 0.1 cm“^ (see 
Supporting Information). 

For the oMX mode, we also find that the surface modes 
have a slightly higher frequency than the modes that are 
localized on the inner layer. Supported by the results 
from our bond-polarizability model (see Supporting In¬ 
formation), we assign the dominant high-frequency oMX 
subfeature to the surface mode, while the faint lower- 
frequency shoulder is assigned to the inner modes. We 
note that the low-frequency shoulder is not observed 
at all at Ei, = 2.33 eV (see Supporting Information) 
and may thus arise from a resonance effect [TOl HH |40] . 
We also point out that ab-initio calculations (see Sup¬ 
porting Information) yield a sizable Davydov splitting 
of 5.8 cm“^ between the bulk A 2 u and i? 2 g frequencies 
(in the absence of the Lydanne-Sachs-Teller interaction). 
In order to reproduce this splitting, an additional force- 
constant between Mo atoms of neighboring layers would 
need to be introduced (which would not modify the split¬ 
ting between the oX modes Bi^ and Aig modes because 
those modes do not involve motion of the Mo atoms). 
Since the resulting fan-diagram is not visible in our ex¬ 
perimental spectra of the oMX mode - due to almost 
vanishing intensity - we did not include this additional 
force constant here. 

Finally, we comment on the possible influence of res¬ 
onance effects on our measurements. Resonance ef¬ 
fects should not impact the phonon frequencies but may 
strongly alter the integrated intensity of one given Raman 
feature and the repartition of the spectral weight within a 
given subfeature. For the iX mode, our bond polarizabil¬ 
ity model predicts spectral weights that are in-line with 
group theory predictions and our experimental results 
(see Supporting Information). However, the same model 
predicts that the highest frequency oX mode should have 
a much larger oscillator strength than its lower frequency 
counterparts. This prediction is clearly in contradiction 
with our observations at = 1.96 eV (see Fig 
where the oX-mode subfeatures have comparable inte¬ 
grated intensities, but is consistent with the absence of 
splittings seen in numerous studies on M 0 S 2 using visible 
photon energies p Ha Eg EDI HD]. Noteworthy, a pro¬ 
nounced splitting of the oX feature was also observed at 
El = 2.33 eV (see Supporting Information). However, in 
this case the high-frequency peak largely dominates the 
oX-mode feature for N > 6. These intriguing observa¬ 
tions, together with the observation of the iX and oMX 
bulk modes, and a particularly intense oX-mode feature 
in monoloayer MoTe 2 at El = 1.96 eV (see Supporting 
Information and Ref. 137!) provide a strong impetus for 
a quantitative analysis of resonant, symmetry-dependent 
exciton-phonon coupling [TDl mi SOI Sni |17j in iV-layer 
MX2. 

Conclusion Using micro-Raman spectroscopy, we 
have reported a unified description of the optical 
phonon modes in a Wlayer 2i7-transition metal dichalco- 


genide crystal (here, MoTe 2 ), between the bulk (three- 
dimensional) and monolayer (quasi-two-dimensional) 
limits. The manifolds of low-frequency interlayer shear 
and breathing modes, and of the mid-frequency modes 
involving out-of-phase intralayer motion of the chalco- 
gen atoms are well understood using classical theories 
of coupled oscillators. In contrast, the behavior of the 
high-frequency modes that involve out-of phase motion 
of the metal and chalcogen planes is largely influenced 
by surface effects. We have introduced a global fitting 
procedure based on a linear chain model to derive the 
force constants up to the second nearest neighbor and to 
assess the strength of the surface effects. This model al¬ 
lows us to deduce the frequency of optical phonons that 
are silent in bulk crystals, namely the low-frequency in¬ 
terlayer breathing mode (LBM) with i? 2 g symmetry, and 
the mid-frequency in-plane (iX) and out-of-plane (oX) 
modes with E 2 U and Biu symmetry, respectively (see Ta- 
ble|!!]). Finally, our work may motivate related studies of 
Davydov splitting, force constants and surface effects in 
other layered crystals, such as black phosphorus [48], bis¬ 
muth selenide and bismuth telluride |49j . 

Methods IV-layer MoTe 2 crystals were prepared by 
mechanical exfoliation of commercially available bulk 
crystals (2D semiconductors) onto Si wafers covered with 
a 90-nm or 285-nm-thick Si02 epilayer. The number 
of layers was first estimated from optical contrast and 
atomic force microscopy measurements, and further char¬ 
acterized using a home-built micro-Raman setup. Micro- 
Raman scattering studies were carried out in ambient 
conditions, in a backscattering geometry using using a 
monochromator equipped with a 2400 grooves/mm holo¬ 
graphic grating, coupled to a two-dimensional liquid ni¬ 
trogen cooled charge-coupled device (CCD) array. Two 
laser photon energies {El = 2.33 eV and El = 1.96 eV 
were employed. Spectral resolutions of 0.6 cm“^ and 
0.4 cm“^ were obtained at El = 2.33 eV, El = 1.96 eV, 
respectively. A laser intensity below 50 kW/cm^ was 
used in order to avoid photoinduced damage of our sam¬ 
ples. In order to attain the low-frequency range (for mea¬ 
surements at El = 2.33 eV), a combination of one narrow 
bandpass hlter and two narrow notch filters (Optigrate) 
was used. After optimization, Raman features at fre¬ 
quencies as low as 4.5 cm“^ could be measured. Polar¬ 
ized Raman studies are performed using a Glan-Thomson 
analyzer placed before the entrance of our spectrome¬ 
ter. An achromatic half-wave plate was placed after the 
analyzer and adequately rotated such that the Raman 
backscattered beam enters the spectrometer with a fixed 
polarization. Thereafter, we will be considering parallel 
(XX) and perpendicular (XY) polarizations of the in¬ 
coming and Raman-scattered photons. The observation 
of residual features from the out of plane modes (LBM, 
oX, oMX) in the XY configuration (see Figs. and. is 
due to depolarization effects induced by our high numer¬ 
ical aperture objective (100 x, NA = 0.9). Finally, the 
measured Raman features are fit to Voigt profiles, taking 
into account our spectral resolution. 
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Supplementary Information 

This document contains the following sections: 

• Force constant model (Section [sTj) 

• Bulk phonon frequencies (Section [S^ 

• Normal mode displacements (Section [S3| 

• Empirical bond polarizability model (Section [S4|) 

• Raman Spectra from the bond polarizability model (Section |S5[ ) 

• Ab-initio bulk phonon frequencies (Section |S 6 [ ) 

• Additional Raman measurements (Section |S7[) 


SI. FORCE CONSTANT MODEL 


As explained in the main text, A-layer MoTe 2 is modeled as a one-dimensional finite linear chain composed of 2N 
Te atoms of mass per unit area and N Mo atoms of mass per unit area (see Fig. [5]. Within one MoTe 2 
layer, nearest neighbor Mo and Te atoms and the pair of second nearest neighbor Te atoms are connected by springs 
with force constants per unit area a and b respectively. Interlayer interactions are described by two force constants 
per unit area /3 and 7 between nearest neighbor Te atoms belonging to adjacent layers and between second nearest 
neighbor Mo and Te atoms, respectively. To account for surface effects, we consider effective force constants Oe and 
(5e for the first and layer. Since substrate effects have been shown to have a negligible influence on the Raman 
modes of MX 2 piiniiiiiiiiisi, we assume that the two extreme layers are only connected to one layer, i.e., we do 
not include an additional spring constant that would account for coupling of one of the outer layers to a substrate. 
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FIG. SI. Schematic of the finite linear chain model. pM (jJ-x) is the mass per unit area of the Mo (Te) atom, a and /? are the 
(intra-plane) force constants that connect the first nearest neighbor atoms. 7 and <5 are the (inter-plane) force constants that 
connect the second nearest neighbor atoms. Uij is the displacement, with respect to the equilibrium position, of the i**' atom 
(i = 1, 3 for Te and i = 2 for Mo) in the MoTe 2 layer (j £ Jl, N|). 


We note Ui^ the displacement, with respect to the equilibrium position, of the atom (i = 1,3 for Te and i = 2 
for Mo) in the MoTe 2 layer {j G |l,iV]). We can then write the equations of motion using Newton’s law. These 
equations form a system of 3N coupled differential equations that can be written as 


dt^ 


-VU, 


(SI) 


M2.1, M3,l, ..., Uij, 


..., 'Ul.JV, U2^N, 


U3,N 


with the displacement vector U = 


and the 3N x 3iV dynamical matrix 
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To find the normal modes, one has to seek for sinusoidal solutions. For this kind of solutions, Eq. (SI I becomes 

VU = u?U. (S2) 


Therefore, the iN normal modes, with eigenfrequencies uJk and normal displacements {k G |1, 3iV]), are obtained 
by diagonalizing the dynamical matrix V. 


S2. BULK PHONON FREQUENCIES 


To obtain the frequencies of the six bulk normal modes, we use the same model as in the previous section [ST] except 
that we apply the Born von Karman periodic boundary conditions to take into account the infinite size of the crystal. 
In this case, the unit cell of this one-dimensional Bravais lattice contains the three atoms of one layer. For the 
layer, we suppose that the equilibrium positions are na for the Mo atom and na — d and na + d for the two Te atoms. 
Thus, Mo atoms belonging to adjacent layers are separated by a. With the same notation as in section we seek 
for solutions in the form of a plane wave with frequency uj and wave vector k : Uj^n = where j = 1, 3 

for Te and j = 2 for Mo, and Aj are constants to be determined, whose ratio specify the relative amplitude and 
phase of vibration of the atoms within each layer. By substituting „ into the equations of motion, we obtain three 
homogeneous equations in terms of Aj. These equations will have a non-zero solution provided that the determinant 
of the coefficients vanishes. This yields 


ika 


[^JLxlo^ - {a + fi+ ^ + 5)Y - 2 (q! -I- 7)] 

-b (a -b 7e®'=“)2((5 -b / 3 e-®'=“) -b (a + -b / 3 e® 

- - 2(a + 7 )] {5 + 

- 2 [a^xw" - (a + /3 + 7 + d)] (a + 7e®'=“)(a + 76-®'®“) = 0 . (S 3 ) 


The Born von Karman boundary condition leads to N nonequivalent values of k given by fc = with p an integer. 
Eq. (S31 does not need be solved for every k. In fact, for the six bulk normal modes, the displacements of the three 


atoms within one layer are either in-phase or out-of-phase with the displacements of the atoms of adjacent layers. 
Therefore, k = 0 oy k = ^ respectively. By Solving Eq. ( |S3[ ) with fc = 0 and k = '^ and using the symmetry of the 
atomic displacements, we can get the expression of the six bulk frequencies associated with the low- (LSM, LBM), 
mid- (iX and oX) and high-frequency (iMX, oMX) modes [ 3 ]. 
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(5 6 ) 

(57) 
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(59) 


Erom the value of the force constants extracted from the fit of our experimental data (see Table II of the main 
text), we notice that |a| » |/3| , I 7 I, |d|. Thus, we can perform Taylor developments of Eqs. (S5) and (S9| to get more 
convenient expressions 
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where fi = 2fj,x + /xm is the mass per unit area of the unit cell. The relative difference between the results of Eqs. 
(S10)/(S11| and the exact values obtained using Eqs. (S5|/(S9|, respectively, is lower than l%o. 

Sri interesting quantity than can be deduced from the expressions of the bulk frequencies, for low-, mid- and high- 
frequency modes, is the bulk Davydov splitting Aw = — uj~ . Again by performing Taylor expansions, we get the 

following expressions for the Davydov splitting 


Awiow 

AWmid 


Awhigh 


P (l_ 7 + 1 ^ 

a V 2a J y ^ix 

f 7^ 

\ /i^ a jj? a J 


a/i 


(512) 

(513) 

(514) 


The deviation of the results of Eqs. (S12) and ( S13 ) from the exact values deduced from Eqs. (S4)-(S7) is lower than 
1%, and the deviation of the results of Eqs. (S14) from the exact values deduced from Eqs. (881 and (891 is lower 
than 10 %. 

Interestingly, the high-frequency Davydov splitting (Eq. (8141) is the only one that can be negative since a >> 7 , <5 

■ 4^ 7 — ^ splitting is normal and the bulk high-frequency in-phase mode has a lower frequency than the bulk 

high-frequency out of phase mode. Otherwise, the splitting is anomalous, as it has been reported for the iMX mode 
in bulk transition metal dichalcogenides [ 20 ] . 


S3. NORMAL MODE DISPLACEMENTS 

Figures [S^ to [S7| show the normal mode displacements associated with the L 8 M, iX, iMX, LBM, oX, and oMX 
modes in A-layer MoTe 2 . 
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FIG. S2. Calculated normal displacements associated with the LSM in = 1 to = 6 layers MoTe 2 . The size of the arrows 


is proportional to the amplitude of of the normal displacement obtained from the solution of Eq. (SI I. The frequencies of 
the modes increase from left to right. The irreductible representation of each normal mode is indicated. The modes that are 
Raman-active in our geometry appear in black. The other modes appear in grey. 
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FIG. S3. Same as Fig. |S2| for the iX modes in A'^ = 1 to A'^ = 6 layers MoTe 2 . 
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FIG. S4. Same as Fig. |S2|for the iMX modes in = 1 to = 6 layers MoTe 2 . 
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FIG. S5. Same as Fig. |S2| for the LBM in A'’ = 1 to = 6 layers MoTe 2 . 
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FIG. S6. Same as Fig. |S2| for the oX modes in = 1 to = 6 layers MoTe 2 . 
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FIG. S7. Same as Fig. |S2| for the oMX modes in = 1 to = 6 layers MoTe 2 . 
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S4. THE EMPIRICAL BOND POLARIZABILITY MODEL 


In the Placzek approximation, the Raman intensity for non resonant first order scattering is given by 

F oc |ej • A'' ■ —(n^ + 1), 

iOy 


(S15) 


where and are the polarizations of the incident and the scattered light respectively, uji, is the Raman frequency, 
and riv = [eyi];>{huj^/U bT) — 1]“^ with T being the temperature. The Raman tensor A^ is related to the change of the 
polarizability a under atomic displacement. It can be calculated as 



(S16) 


where is the Cartesian component of atom 7 of the orthonormal vibrational eigenvector My is the atomic 
mass, and 


dEidEjduky dukj ' 


(S17) 


E is the total energy in the unit cell, ill is a uniform electric field, Uk-y is the fc**' component of the atomic displacement 
u of the atom 7 and is the electronic polarizability tensor. 

The bond polarizability model approximates the total polarizability of the atom as the sum of the individual bond 
polarizabilities. Under the assumption that the bonds have cylindrical symmetry, the polarizability tensor for a single 
bond can be written as 


^ij — ^ T E {cXl ^p) ( RiRj o 

CXpSj^j E (cT/ CXp^R^Rj 


(S18) 


where ai and ap are longitudinal and perpendicular polarizabilities of the bond and Ri are the components of the 
unit vector along the bond. Moreover, in the model, the polarizabilities depend only on the length of the bond 
{R = VR^). The bond vector R joining atoms 7 to atom 7 ' is given by 


and 


Rj - Ry Ry' ‘Uy 'Uiyf 


(S19) 


dajR) 

duky 


da{R) dR 
dR duk-y 
da{R) 1 

dR 2^/^UM^ 

dR 2R^ ’ 


R 


9R 

duk^ 


5R 

duky 



(S20) 


Hence, 


da{R) 

dUky {1) 


= Cx'Rk 


where a' =- ^ . The contribution of a particular bond b to the B tensor is therefore 




du 


k-y 


with 


— a'RkSij E {ct'i — a')RiRjRk E {cxi — ap){{dkRi)Rj E {dkRj)Ri) 


dkRi — {^ik RiRk 


(S21) 


(S22) 


(S23) 


In the 2E[ structure of MoTe 2 , the Molybdenum atom is bonded to six Tellurium atoms and the Tellurium atom is 
bonded to three Molybdenum atoms. In our calculations, we neglect the weak inter-layer bonds. That is the reason 
why we don’t get a finite peak for the E^g (shear) mode with this model. 

The calculated B tensors are as follows: 
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The values of the constants a, b, p and q, in terms of the polarizabilities and their derivatives, obtained after substi¬ 
tuting the MoTe 2 lattice parameters are: 

a = —0.24q;; — 0.56a; -I- 0.24a; — 1.40a; 
b = 0.48a; - 0.84a; - 0.48a; - 1.12a; 
p = 0.28a; -I- 0.64a; — 0.28a; — 0.64a; 
q = 0.18a; — 0.56a; — 0.18a; -I- 0.56a; 


(S24) 
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S5. RAMAN SPECTRA FROM THE BOND POLARIZABILITY MODEL 


With the bond polarizability model, we can assign a Raman intensity to each Raman frequency obtained with the 
force constant model. While active Raman modes agree with group theory and experimental data, the model does not 
include Raman resonance effects, and some discrepancies arise when we compare with experimental spectra. Figure 
S8 shows the theoretical Raman spectra for the modes iX, oX, iMX, and oMX. The spectra for different n values are 
offset for clarity. The iX mode has \N/2\ active subfeatures in A—layer MoTe 2 and is in nice agreement with the 
spectra of Fig. 3. 

For the iMX mode, there are two close peaks for A > 3; one being an inner mode (with lower frequency) and the other 
being a surface mode (with higher frequency). The difference in their frequencies is around 0.5 cm“^.The absolute 
intensity of the inner mode increases almost linearly with number of layers, as there are more layers vibrating, whereas 
the absolute intensity of the surface mode is independent of number of layers, as only the outer molecules are vibrating. 
As a result, the relative intensity of the surface mode drops as A increases, and thus the maximum of the combined 
peak shifts to smaller frequencies as A increases. 

In the case of the oX mode our model reproduces the [A/2] active subfeatures well, but it fails in describing the 
observed relative intensity between peaks. The highest frequency peak in the oX-mode feature has much larger relative 
intensity than the rest of the peaks in the model, which is in contrast to our measurements, where all the phonons 
have comparable intensities. We can understand this result from the shape of the tensors shown above. The 
largest contribution occurs when all the Te layers are out-of-phase with the neighboring Te layers. In case two layers 
vibrate in-phase, they cancel each other, and therefore increasing the number of layers in-phase reduces the intensity 
drastically. The similar intensity of the A] peaks suggest the participation of Raman resonance effects. 

For the oMX mode, the surface phonons, which are split from the inner modes, have much larger intensity than the 
inner modes in the model. From the Raman tensors, if tellurium atoms in the same layer vibrate in-phase and with 
identical amplitude, the intensity is identically zero. Only differences in the amplitudes can generate Raman signal, 
which will be very small in comparison with the intensity of the oX modes. The largest amplitude differences are 
due to surface phonons, which carry most of the contribution to the Raman intensity and it results in the sole peak 
observed in the Raman spectra. A feature assigned to the inner modes is slightly visible in the measured spectra at 
El = 1.96 eV (see Fig. 3(c) in the main text) but not at El = 2.33 eV, which, again, suggests a resonance effect at 
El = 1.96 eV, similar to the oX mode. 
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FIG. S8. Raman Spectra obtained with the bond polarizability model. 
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S6. AB-INITIO BULK PHONON MODES 


We have complemented the empirical model with ab-initio calculations of the phonon modes for bulk MoTe 2 . We 
have used density functional perturbation theory (DFPT) [50] as implemented in the Quantum Espresso code [ST] , 
We used the local-density approximation (LDA) which does not properly take into account van der Waals interaction 
between the layers but nevertheless gives decent result for the phonons of many layered systems because it over¬ 
estimates the weak covalent part of the inter-layer bonding. The energy cutoff is 80 Ry, and the Monkhorst-Pack 
sampling of the k— grid is 12 x 12 x 4. The optimized lattice vectors are a = 3.499 Aand c = 13.829 A. Figure S9 
shows the ab-initio frequencies together with the corresponging phonon eigenvectors. We have grouped the phonon 
modes in Davydov pairs, except the shear and layer-breathing modes (whose “Davydov partner” would be a zero 
frequency acoustic mode). We observe a positive Davydov splitting for the iX, oX and oMX modes and a negative 
Davydov splitting for the mode iMX mode, as reported for other transition metal dichalcogenides [5D]. The calculated 
Davydov-splitting of the iX mode is 2.1 cm“^ (compared to the 2.7 cm“^ extrapolated from the measurements). For 
the oX mode, we obtain a splitting of 3.2 cm“^ (compared to 4.7 cm“^ from the experiments). The agreement between 
calculations and experiment is not perfect because of the lack of a proper treatment of van der Waals interaction 
with local exchange-correlation functionals. The predicted Davydov-Splitting for the oMX mode is 5.8 cm“^. For the 
calculation of the frequency of the A 2 u mode, we have not taken into account the coupling to an external-electric field 
(Lydanne-Sachs-Teller interaction) that would lead to an up-shift of this mode. However, in layers of finite width, 
this interaction is absent. 
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FIG. S9. Optical phonon modes of bulk MoTe 2 . In the first row, modes with in-plane polarization in ascending order of 
frequency. In the second row, the out-of-plane modes are shown. Davydov pairs of phonon modes are plotted in one box. 
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ST. ADDITIONAL RAMAN MEASUREMENTS 


Figure [S10| shows the raw Raman spectra of A^-layer MoTe 2 recorded at i?L = 2.33 eV and = 
that the iX mode has not been studied at E]^ = 1.96 eV due to the relatively large bandwidth of our 
Ei^ = 1.96 eV. 


1.96 eV. Note 
Notch filter at 


(a) 


= 2.33 eV 



(b) El = 1.96 eV 



FIG. SIO. Micro-Raman spectra of A—layer MoTe 2 recorded under the same conditions at a photon energy of (a) 2.33 eV and 
(b) 1.96 eV. The spectra are vertically offset for clarity. 
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Figure Sll shows Raman spectra of the oX, iMX and oMX modes in TV-layer MoTe 2 recorded at Ei, = 2.33 eV. 
The results recorded at i?L = 1-96 eV are discussed in the main manuscript. At Al = 2.33 eV, the Davydov splitting 
also appears clearly for the oX feature, although the highest energy subfeature contains most of the oscillator strength 
for N > 6. The iMX feature also downshifts as N increases and no appreciable splitting can be resolved. However, in 
contrast with our results at Al = 1-96 eV, the oMX feature does not exhibit any measurable splitting at Al = 2.33 eV. 
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FIG. Sll. Normalized Micro-Raman spectra of the (a) oX, (b) iMX, and (c) oMX mode-featnres in A—layer MoTe 2 recorded 
at Al = 2.33 eV. The measured Raman features (symbols) are fit to Voigt profiles (solid lines). For the modes that show a 
Davydov splitting, each subpeak is represented with a colored dashed line. A featureless background (grey dashed line) has 
been considered when necessary. 
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The corresponding fan diagrams associated with oX-, iMX and oMX-mode frequencies recorded at i?L = 2.33 eV 
are shown in Fig. |S12[ together with the fan diagrams for oX and oMX modes extracted from the data recorded 
at i^L = 1-96 eV and discussed in details in the main text. These two sets of data are very consistent with each 
other. Still, we can notice a small rigid shift of approximately 0.2 cm“^ which is smaller than the resolution of our 
experimental setup. This shift presumably arises from uncertainties (below our spectral resolution) in the calibration of 
our spectrometer. Importantly such a small shift has a negligible influence on the determination of the force constants. 
Indeed, the latter vary by less 1.5 % if one uses the oX, oMX and iMX frequencies recorded at Ei^ = 2.33 eV instead 
of their values recorded at Ei, = 1.96 eV in the global fitting procedure described in the main manuscript. 


oX iMX oMX 



N N N 


FIG. S12. Frequencies of the (a) oX, (b) iMX and (c) oMX modes extracted from Voigt fits as a function of the number of 
layers N. Green squares (red circles) correspond to data recorded at E-l = 2.33 eV {E-l = 1.96 eV). 



















